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We develop a framework for constructing initial data sets for perturbations about spherically 
symmetric matter distributions. This framework facilitates setting initial data representing astro- 
physical sources of gravitational radiation involving relativistic stars. The procedure is based on 
Lichnerowicz- York's conformal approach to solve the constraints in Einstein's equations. The cor- 



^ • respondence of these initial data sets in terms of the standard gauge perturbation variables in the 

I Regge- Wheeler perturbation variables is established, and examples of initial data sets of merging 

. neutron stars under the close-limit approximation are presented. 
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I The early years of the next millennium will hopefully be remembered for the birth of Gravitational- wave Astronomy. 
\ With several large scale interferometers (LIGO, VIRGO, GEO600, TAMA) under construction, and the continued 
0^ improvement of the technology for cryogenic resonant-mass detectors (ALLEGRO, AURIGA, EXPLORER, NAU- 
0]S TILUS), there are many reasons to be optimistic at the present time. However, the interpretation of data from the 
^ ^ new generation of detectors will heavily depend on accurate "templates" of gravitational wave-forms. For a given 
astrophysical source of gravitational radiation, construction of such templates involves fully non-linear or perturbative 
I \ approximations to Einstein's field equations. In both instances, the construction of appropriate initial data consti- 
bJO' tutes a fundamental issue. It is absolutely necessary that initial data represent a "realistic" stage of the astrophysical 
^ I system under consideration. 
• rH , The early years of Numerical Relativity were in part characterized by studies aimed at constructing initial data for 
Einstein's equations, namely data that satisfies the Hamiltonian and momentum constraints. Of particular interest 
was obtaining solutions to the constraints which represent black hole binaries [Q . These initial data studies highlighted 
the importance that Lichnerowicz- York's conformal approach |^ plays in facilitating solving the constraints. 

In the perturbative arena, initial data sets also play a fundamental role if a connection with systems of astrophysical 
relevance is to be made. Examples of perturbative studies where astrophysically consistent initial data is needed are 
the point-particle [|| and close-limit approximations [Q to black hole coalescences. For the close-limit approximation 
in particular, this issue is crucial. The focus is then on the late stage of the merger, when the binary system can be 
approximated as a single, perturbed black hole. In the case of perturbations of relativistic stars, most of the studies 
have not considered initial data with direct connection to a given astrophysical situation. In other words, the focus 
has so far been on investigating how the star reacts to generic perturbations. 

The goal of this paper is to provide a mechanism for generating initial data for perturbations of a relativistic star. 
This data should ideally represent astrophysical situations of relevance to gravitational- wave detectors. Inspired by 
the success that Lichnerowicz- York's conformal approach has enjoyed for solving Einstein's constraints, we base our 
methodology on "linearizing Lichnerowicz- York's procedure" . By doing so, we take advantage of the prescription for 
knowing which pieces of information, among the metric and its "velocity" , are fixed by the constraints and which 
are freely specifiable. Moreover, we inherit the machinery used in the past for the construction of initial data sets 
representing binary systems. That is, the procedure described in this paper provides a natural framework for obtaining 
initial data in connection with the close- limit approximation to neutron star mergers 1^. As with black hole binaries, 
the close-limit to neutron star collisions deals with the late stages of the merger, at the point in which the systems 
can be approximated as a single neutron star "dressed" with perturbations. 
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II. LINEARIZATION OF LICHNEROWICZ- YORK'S CONFORMAL APPROACH 



Given a spacetime with 4-nietric g^i, and a foliation of this spacetime with Eulcrian observers having a 4-velocity 
fi'^, the constraints in Einstein's fields equations can be written in a 3+1 (or ADM) form Q as 

R + K^ - k,jk''^ = 16 TT CT (Hamiltonian) (1) 
Vj{k'^ --h'^k) =8TTf. (Momentum) (2) 

Above, hij is the 3-metric (/ip,y = 5^,^ + h^h^) and Kij the extrinsic curvature of the time-like hypersurfaces in the 
foliation, with K — h^^ Kij and R the scalar curvature. Furthermore, is the covariant derivative associated with 
the 3-metric hij , and a and are the energy and momentum densities of the matter sources. Greek letters denote 
spacetime indexes and Latin letters spatial indices, and we use units in which G = c = 1. Contrary to the common 
practice, we use "hats" to denote physical space since we will be mostly working in the conformal space. 
For a perfect fluid, the stress-energy tensor is given by 

Tfj^L' ^ {p + p)u^_,Uu+p g^ij , (3) 

where p and p are respectively the total mass-energy density and pressure of the fluid measure by an observer with 
4-velocity ii^. The energy and momentum densities appearing in the constraints are obtained from 

a ^f^.nf'n'' = (p + p) f ^ p (4) 
f = -f^^n''h'''^'={p + p)Jh^"'u^ (5) 

where 7 = is the relativistic boost factor. 

The fundamental virtue of Lichnerowicz- York's conformal approach |^ for solving the Hamiltonian and momentum 
constraints is that it provides a concrete recipe for singling out which four "pieces" among the twelve components 
{hij, kij) are to be solved from Eqs. (|l|) and (||). The starting point is the Ansatz 

hij = (t>'^hij, (6) 

where the conformal metric hij is assumed to be known. Thus, for the metric, the piece that is fixed by the constraints 
(Hamiltonian) is the conformal factor cj). The other three quantities fixed by the constraints (momentum) involve 
the extrinsic curvature. The idea here is to decompose the extrinsic curvature into its trace, tracefree-transverse and 
tracefree-longitudinal parts. To achieve this, the extrinsic curvature is first split into 

K'j ^ A'j + l^h'^k. (7) 

Before decomposing the tracefree part A^^ into its transverse and longitudinal parts, the following conformal trans- 
formation is applied: 

1^=0-10.4*^'. (8) 

The exponent in the conformal transformation is motivated by the fact that this transformation possesses the following 
property: 

VjA'^ ^(t)-^^VjA'^, (9) 

with Vi covariant differentiation associated with the background metric hij. This property simplifies the conformal 
transformation of the divergence of the extrinsic curvature in the momentum constraint (|^) . 

Since, as will become clear below, the trace of the extrinsic curvature K is not fixed by the constraints, no conformal 
transformation is imposed on the trace of the extrinsic curvature {K = K) . Once the conformal transformation is 
applied, the next step is to decompose Aij into its transverse and longitudinal parts, namely 

A'^ = A^^ + {lWf\ (10) 

where 
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{IWy^ = 2\/^W^ - ^/l*iVfcH/^ (11) 

v,a:^' - 0. (12) 

With the above conformal transformations and transverse-longitudinal decompositions, the Hamiltonian and momen- 
tum constraints become 

8V*V,;(/.- + - ^K^<j)'^ + 16 ncrcf)-^ = (13) 

y^VjW + ^V'VjW^ + R'jW^ ~ V'K - 8 TT = , (14) 

where is the 3-Ricci tensor of the conformal space and R its trace. In deriving Eqs. (|l^ ) and ([T^), the following 
conformal transformations for the energy and momentum densities were used: 

a = <ir^a (15) 
f^r^^f. (16) 

To summarize, the Hamiltonian constraint fixes the conformal factor and the momentum constraint determines 
the generator VF* of the longitudinal part of the conformal, traceless part of the extrinsic curvature. The freely 
specifiable data in this coupled set of equations are the conformal metric /i^, the trace of the extrinsic curvature if, 
the source functions (ct, j'*), and the divergence free, traceless part of the extrinsic curvature A^-*, which is hidden in 
Ai, in Eq. (|i|). 

Thus far, we have just reviewed Lichnerowicz- York's treatment of the initial data problem. We now introduce our 
jirst assumption. The initial data, (ft-ij, A*-', if, o", f \ are assumed to be close to a given background, i.e., 

(17) 
(18) 

k = if (0) + if(i) (19) 

0- = CT(o) + <T(1) (20) 

/=i|o)+ia)' (21) 

where labels background quantities and ^^-^ first-order perturbations. The corresponding decomposition of the 
conformal space quantities yields: 

(22) 

(i, = + (23) 
J^i^J^l^^^l^ (24) 
= + (25) 
if = if + if (26) 
a = (27) 

(28) 



where 





/ ~ J(0) +■?'(!)' 




ii(°) 


- <P(0) 






= -^(0) hf^+^4>%) hi) 




^(0) 






^(1) 




^(1) ^(0) 








^(1) 


=^(1) 





(29) 
(30) 
(31) 
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o-(o) = 


^(0) ^(0) 




o-(i) = 




CT(0) 




1-10 -i 
^(0) J(o) 




^(1) = 




(1) J'co) 



(35) 
(36) 
(37) 
(38) 

Although aU scalars, vectors and tensors are expanded in order of smaUness, it is important to stress again that all 
but (j) and are freely specifiable, and this property is independent of the order of the perturbation. 

At this point we introduce our second assumption, which is that the perturbations of the conformal background 
vanish on the initial data slice; that is, 0. The primary motivation for this choice is the simplification of the 

coupled system of constraint equations. We shall later discuss the physical relevance, as well as the implications and 
restrictions, that this assumption imposes on the class of initial data that one has access to with our procedure. Using 
h\^^ — 0, together with the above perturbative expansions in the coupled elliptic system (^ and (p^), we obtain 



SV^V 



i<p(l) 



R + TA^^^Af ct>l,] + yi^fo) <I>U + 48^a(o) ^^^^ 



(1) ^•'•'(1) gV'CO) 



«(o) '^'(i) 



'^J(i) 



(39) 
(40) 



where R, Rij and refer to the background. In writing Eqs. (|3^) and (^), we have used the fact that the zeroth-order 
quantities satisfy the constraints. Notice that, as in the non-linear case, the constraints remain coupled. 



III. THE INITIAL DATA PROBLEM FOR RELATIVISTIC STELLAR PERTURBATIONS 

We now focus on constructing initial data sets for which the background is a static and spherically symmetric stellar 
model, with 4-metric given by 

ds^ = -e^" dt^ + e^^ df^ + dO^ + sm^ 6 dip^ , (41) 

where the metric coefficients v and A are functions of the radial coordinate f only. Einstein's equations for this 
background reduce to solving three equations. The first equation defines the "mass inside radius f; 

^=4^rV(o). (42) 
Here m is a function of f that is related to the metric function A by 

T\ 2 TO , , 

e-^^ = l —. (43) 

r 

Equation (|4^) is directly obtained from the Hamiltonian constraint. The second equation belongs to Einstein's 
evolution equations and reads 

dv 

^ = 7^("^ + 4^^'P(0))- (44) 
Finally, conservation of momentum yields the condition for hydrostatic equilibrium: 

^ = -(p(0)+P(0))^ (45) 

The above stellar structure system of equations (the Tolman-Oppenheimer-Volkoff equations) must be supplemented 
with an equation of state. For simplicity, we use the polytropic equation of state p — np^, where k and F are the 
adiabatic constant and index, respectively. The adiabatic index F is related to the polytropic index n by F = 1 + l/n. 
In the specific example provided later, we use F = 2 (n = 1). 
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We now assume that the zero-order quantities (/5(o), ^(0), A, v) have been obtained from solving Eqs. (^2[), (H) and 
(^). The next step is to solve for the first-order perturbations from the linearized constraints (|39| ) and ([40|). To 
facilitate this task, it is convenient to perform a coordinate transformation and bring the 3-metric in (^Tj) into the 
isotropic, conformally flat form: 

ds^ 0^0) (^^^^ + + sin^ 6* ) , (46) 

where the conformal factor is given by 

<^(0) - [-) (47) 

and the transformation of the radial-coordinate is obtained from 

dr ^ r , 

(48) 

dr r 

The static nature of the background spacetime and the gauge choice of a vanishing shift vector imply if ^ °^ — j^^^^ ~ 
for this background spacetime. This leads to considerable simplifications in the following, but the procedure we 
describe for constructing perturbative initial data can easily be extended to include also time-dependent, spherically 
symmetric background spacetimes. 

We introduce here our third assumption, which is that ^*"(]^) = ^(i) — 0- The vanishing of the transverse-traceless 
part of the perturbation to the extrinsic curvature has direct implications to the gravitational radiation content of 
the initial data. In a way this assumption can be viewed as "minimizing" the amount of gravitational waves in the 
initial spacetime. That this may not be desirable is obvious, but in order to be able to assign this part of the free data 
physically correct values we need a detailed knowledge of the past history of the system. Such information requires long 
term, nonlinear evolutions and is far beyond our present capabilities. The present assumption is convenient in that 
it simplifies the calculations considerably. Furthermore, we are unlikely to overestimate the amount of gravitational 
radiation emerging from true physical systems if we base our estimates on the present approach. The vanishing of the 
trace of the extrinsic curvature to first order is less restrictive physically. It simply implies (using also the fact that 
^(0) = 0) that the slicing of the perturbed spacetime is maximal to first order. With those further assumptions, the 
linearized constraints (39) and ( ^ ) decouple and take the form 

V-'Vj(/)(i) = 6 TT (7(0) 0(o') ~ 2 7r(?!)jr^^cr(i) (49) 
W^U = W.,V' (50) 
V^V,F^ = 8^j(i). (51) 



where we have decomposed the vector W^(\') in Eq. (^) foUowini 

1 
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^V- -V'U. (52) 



From now on, we will drop the label (i) in j^'j^^ and since the zero-order values for these quantities vanish. 

The set of equations (|49|-^l[) constitute a coupled set of elliptic equations in three dimensions, expressed e.g., in 
the coordinates (r, 6, ip) of the background space. Due to the spherically symmetric nature of the background, these 
linearized equations allow for a separation of variables. Specifically, we can apply a spherical harmonic decomposition 
of form: 

(53) 
(54) 
(55) 
(56) 
(57) 
(58) 







= ^0(i)(r) Yi„,{0,ip) 
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hn 
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= E^i(^)^i(^'^)+''-^2(r)e*2(6 




-T.h{r)e\( 
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where 



0, 0) 

0, \ dgYir, 



sin^ I 



if ^ Im 



(59) 



= 0,- 



1 



■ dgVir. 



Here e\ and 62 are the basis vectors of even-parity perturbations and e\ is the basis for odd-parity perturbations. In 
the above expressions and what foUows, it is understood that the radial functions are for a given (I, to). These indices 
have been suppressed for economy in notation. We note that in terms of the radial functions, Eq. ( p2[) reduces to 



Wi=Vi^]-^U 
4 dr 

W2 = V2--U 
4 

w. = v.. 



(60) 



After separation of variables, the system of equations (^9|-^ is rewritten as a system of coupled radial elliptic 
equations: 



1 d 

j,2 flj. 



dr 



6n a 



1 d 

j,2 flj. 



dr ( dr^ 
W + 1) + 2] 



(0) 'P(O) 



= -27r 



U 



Vi 



V2 



'^(0) ^(1) 

1 d / 2 y\ 

Stt Ji 



/(/ + 1) 



V2 



l±,r^-V2 
dr \ dr 

]_d_ 

J.2 fj^j. 



2 d . , 
r —V3 
dr 



l{l 



TT J 2 
TT J3. 



(61) 
(62) 
(63) 
(64) 
(65) 



Equations (|6l[-p5D fully characterize, for each (Z,to) harmonic, initial data {(t>(i),U,Vi) to first perturbative order, 
once the background conformal factor 0(q) and density ct^q) as well as the fluid perturbations cr(i) and Ji are specified. 
Outside the sources, the solutions to the above equations are U = ur"", Vi ~ vi r**, V2 = V2 r'^ and V3 — V3 r"* with u, 
vi, V2 and V3 constants. In order to have regular solutions for r — > 00, one needs a = d = —{I + 1) , b = c = —I and 
vi = lv2, where we assume that Z > 1. As we shall later see, the corresponding interior solutions for head-on and 
inspiral close- limit collisions exhibit the same scaling; that is, Vi — IV2 as well as Ji = I J2- With this assumption, 
the system of equations (^-65) reduces to 



c?r i dr ' 



dr \ dr 



j,2 flj. 



dr 



r^ 



dr 



(Z-1) 



Vi 



Vs^SirJs. 



(66) 
(67) 
(68) 



IV. CORRESPONDENCE WITH REGGE- WHEELER VARIABLES 



Before we proceed to present examples of initial data sets constructed using the above approach, we want to 
establish the correspondence between our variables and the standard Regge- Wheeler variables. This is relevant since 
the Regge- Wheeler notation (and the associated gauge) is customarily used in perturbative evolutions for spherical 
stellar models Q. To this end, consider a spatial tensor Tij such that Tij{f, 9, (p) — T^'^\f) -\- T'>j\f, 0, (p). The 
perturbations of this tensor can be decomposed as 
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TS^=ti{r) 4 +H2{f) fl + (r) + U (r) ^ + h (r) + h (f) /f^- , (69) 

where 

/Yi„^ 0\ 

4= 0, (70) 

\ 0/ 



/ dgYhn d^Yi„r\ 

/2. ^ symm , (71) 

\ symm / 



/O 

Vo 



^Y 



hn > 



(72) 



/O \ 

4 = a^Yi^ (9^9^- cot ea^)^;™ (73) 
V symm ((9^ + sin 6* cos Gde)Yi,n) 



4 = symm ) , (74) 

\ symm 




- ii^ - cot e) Yi^ -j^{^^dl + cot 9de - df) Yi^ . (75) 
\ symm — su\9{ded^ — cot 9d^) Yim I 

Above, 4, 4, /? , and 4 represent the even-parity tensor spherical harmonics and 4 ^^^^ 4 the odd-parity 
counterparts. 

Using the Regge- Wheeler notation, the perturbations of the spatial metric (in physical space) read 



^ = E ^2 (r) 4 + /.r"" C'^) 4 + A^f) 4 + f2 G(f) 4 + (f) 4 + (r) 4 , (76) 



where K must not be confused with the trace of the extrinsic curvature. From the previous section, we have that the 
spatial metric can be constructed from 



= h^^ + }^ = <f = (0(0) + 0(1))' 

= c^l,)hf+i4,l,^ct>ii)hf (77) 
= ^!^'+4^/.g), (78) 

where 0(o) — (f/r)^/^, h^^^ — diag(e2'^, f^, sin^ 0) and 0(i) — X]/m 0(i) ('') '^ith 0(i)(r') a solution of the radial 
equation (pl|). Thus, our approach to construct initial data yields spatial metric perturbations of the form 

Im ^ ^ 

Comparison of the metric perturbations ( |79| ) with (|7^) shows that our procedure for constructing initial data yields 
in terms of the Regge- Wheeler notation h^''' = h™'^"' = ft.2 = G = and 

-\ 1/2 



F2 = i^ = 4 0(i)(r) - . (80) 
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Consider now the extrinsic curvature. As mentioned before, the extrinsic curvature vanishes to zero-order. In 
addition, we made the choice of having vanishing first-order trace and transverse-traceless parts Thus, 

the extrinsic curvature is completely determined by the vector and the conformal factor of the background space 
0(0) from 



In terms of tensor spherical harmonics, the extrinsic curvature reads 

= E (^) + ^ ^2(r) 4 -t- r2 fc3(r) ^ + k,{r) ^ -|- kr-,{T) ff^ + k,{r) ffj 

Im 



(81) 



(82) 



where 



fci 

fc2 

ki 



2 

37 
1 
r 
2 

r 



2 r — Il^i + 1(1 + l)W2-2Wi 
ar 



r—W2 + Wi-W2 
dr 



~r— W1 + W1+ 1(1 + 1) W2 
ar 



(83) 



fcs 

ke = -2rW3. 



r W3-W3 
dr 



(84) 



Before proceeding, it is appropriate to discuss the implications of the assumptions we imposed for the presence of 
even and odd parity perturbations in the initial data. 

Consider first the perturbations of the 3-metric: The assumption of vanishing perturbations of the conformal 
3-metric implies (as seen from Eq. ([79|)), that h°'^'^ = hf^" = /12 = G = 0. This means that all odd parity 
perturbations of the 3-metric must vanish on the initial surface. This does not however exclude the presence of odd- 
parity perturbations. Such perturbations may enter via the Regge- Wheeler variable ho, which is freely specifiable. 
Similarly for even-parity perturbations, the Regge- Wheeler quantities Hq and Hi are not part of the constraints and 
can be chosen freely. A choice of these three variables correspond to choosing a slicing for the spacctime, i.e. specifying 
the lapse and the shift vector. Specifically, Regge- Wheeler gauge corresponds to hf"^ = G = /12 = above, as well 
as perturbed lapse 



(85) 



and shift vector 



(86) 



Furthermore, there are no constraints on the odd or even parity nature of the extrinsic curvature initial data. This 
follows immediately from Eq. (|8^ ) since the coefficients ki through fcg are in general non- vanishing. As a consequence, 
even if the initial 3-metric has vanishing odd-parity perturbations, their time evolution will in general include such 
perturbations. 



V. SAMPLE INITIAL DATA SETS: COLLIDING NEUTRON STARS 



We will now apply the method for constructing initial data, that was presented in Sec. Ill, to a case of astrophysical 
relevance. We consider collisions of neutron stars under the close-limit approximation. Related results have already 
been discussed by Allen et al That study was specialized to the case of head-on collisions. Here, we present 
a general discussion of neutron-star close-limit initial data, and give explicit results for both boosted head-on and 
inspiralling collisions. 
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As stated at the end of Sec. Ill, initial data are obtained by solving the system of equations equations (^^^5|). 
The input for these equations are the background conformal factor ^(q) and density cr(o) and, in addition, the fluid 
perturbations (J(i) and Ji. For close- limit collisions, one specifies the background from the outcome of the collisions. 
The fluid perturbations cr(i) and J^, on the other hand, are obtained by "subtracting" the background from a suitable 
superposition of stars that represent the initial configuration. A suitable way of relating the two initial stars to the 
final configuration was presented in and we refer the reader to that paper for further details. 

Let us first consider the perturbation to the background density. Neglecting complicating factors, such as the effects 
from tidal deformations, we approximate the total density of the binary system with the following superposition of 
density profiles of isolated neutron stars Q: 



1/2 



(87) 



where cr* is the conformally transformed density profile of the colliding neutron stars in isolation located a distant 
^* in conformal space. This functional form (|8^) for the total density is chosen since it leads to the correct zero- 
separation and infinite-separation limits. That is, for zero-separation ct — > ct*, and for large separations cr(r*) — > 
(7*(r* — O + cr*(?'* + 'f*). We now introduce the close- limit approximation, ^' <C 1, and write 



Therefore, Eq. (^) takes the form 



(89) 



Given this total energy, the density perturbation is obtained from 

(7(1) = cr - (7(0) 

= ^T* - fT(o) + \ (f V, a^f + i f f V, V, a, . (90) 

Similarly, the momentum density is given by the superposition of momentum densities of boosted isolated stars: 

nr^)^3\{r'' +e)+3-{r'' -e)- (91) 
In this case, there is no need for a counterpart of the last term in (^). For both, head-on and inspiral collisions, 

j+ = -J- =il; thus. 



/(r-')=j:(r'^+e')-j:(r'=-e') 



(92) 



which obviously has the appropriate zero-separation limit, namely i^{r^) — > as f — > 0. Once again, we apply the 
close-limit condition and approximate 



(93) 



Therefore, 



Notice that Eq. 
static. 



f^^e^jji- (94) 

directly gives the momentum density perturbation because, by construction, the background is 



A. Head-on Collision Initial Data 



For simplicity, we assume that the collision takes place along the z-axis. Therefore, 

f — - / 



(95) 
(96) 
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where 



CO8 0, 



sin 9, 



(97) 



is a unit vector along the z-axis. 
Substitution of iM) into Q yields 



O'(l) = cr* - cr(o) + 2 ^ 



cos 



dr 



cos 



io 



6^ 1^;^"* 



dr 



dr 

1 d 

— 

r dr 



■ 2n^ d 
cr, + sm - -r^ 
r dr 



2 d 
r dr 



(98) 



It is clear from (^8[) that the conformal density perturbation has a monopole contribution {m — Q, I — Q). Since our 
main interest is associated with the gravitational waves generated during the merger, we will ignore this contribution; 
it obviously does not lead to any gravitational waves. The remaining part is the radiative quadrupole (m = 0, / = 2) 
perturbation: 



o-(i) 



15 



dr^ 



dr 



1 d 

r dr 



(99) 



This initial data set is discussed in considerable detail in 

Let us now consider the momentum density perturbation. Substitution of ( p6| ) into 

= -2(-^J^ ( cos2 0,-isin6' cos0,O 
dr \ r 



yields 



(^00,0,0) 

3 dr 



15 



dr 



Y2o,-deY2o,0 
r 



(100) 



As expected and in agreement with the density perturbation, we agai n have a monopole momentum density pertur- 
bation that we shall ignore. Comparing the quadrupole terms in (|10C|) with (pq), we have J3 = and 



(101) 



Ji = 2J2 = -— V57r^— J*. 

15 dr 

Fr om E g. (Q) and (^ and the conformal transformations (|l5|) and (|l^), the expressions for ct* and J, in Eqs. (|& 
and (101) are given by: 



J* = (p* -|-j5») u0(o) 



(102) 
(103) 



with V the magnitude of the collision velocity. In writing the above expressions, we used that to zero-order — 
and 7 1 + 0(iL^). Notice that the conformal factor in connection with the background star, (^(0), was the one used 
in Eqs. ( |102| ) and (103) to transform the physical TOV solutions of the colliding stars since it is the background star 
that provides the conformal space where our calculations are performed. Also important is to notice that, once the 
background and colliding neutron stars models have been completely determined, there are only two parameters that 
characterize the initial data: the separation ^ and the velocity v. 



B. Inspiral Collision Initial Data 



For this case, we assume that the initial configuration is such that the neutron stars are along the x-axis and their 
momentum pointing along the y-axis. Then 
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where 



sm U cos ip, - cos a cos Lp, 
r 



sm ip 
r sin0 

cos Lp 



sm C/ sm 93, - cos y sm (/3, 

r r sm 



(104) 
(105) 



(106) 
(107) 



are unit vectors along the x-axis and y-axis, respectively. Substitution of (104) and (105) into (M) yields 



CT(i) 



= (T* - (7(0) + 2 ^^'^ ^ '^^^^ 

+ (cos^ 6' cos^ ip + s\v? p)--^a^ 

r dr 



dr 



dr'^ 



00 



a* - (7(0) + - e ^ + ( - 



dr 



2 d 
7' dr 



Rer22 r 



dr2 



dr 



1 d 
r dr 



(108) 



As with the head-on collisions, we concentrate on the radiative quadrupole {m = 2,1 — 2) term, so the density 
perturbation cr^i) source of Eq. ( |6l] ) is given by: 



^(1) 



dr^ 



dr 



1 d 
r dr 



(109) 



Notice that the only difference between (109) and the corresponding source term in the head-on collision case, i.e. 

Eq. (p^), is a numerical factor. 

For the momentum density perturbation, substitution of (105) into (Q) yields 

f 



2ix' Vji.hy') 
dr \ 



1 



1 



sm p cos p, — sm t/ cos u sm p cos (p, — cos p 
r r 



0,0, 



1 



1 



M^i^J*lm [¥22, -de Y22 , — ^ 
V 15 dr \ r r sm ( 



■d^Y22 



Comparing ( |110| ) with (|58|), we deduce for the quadrupole term (m = 2, ^ = 2) J3 = and 



(110) 



(111) 



The dipole term (m = 0, / = 1) does not contribute to the emerging gravitational radiation and can be ignored. Once 
again, the momentum density perturbation for the inspiral case only differs from the head-on case by a numerical 
factor. The only non-trivial differences in the initial data will then arise fr om t he Yi ^s s ince in one case m = 
(head-on) and for the other m = 2 (inspi ral) . The quantities ct* and J* in ( |l09| ) and ( 111 ) are obtained as in the 
head-on collision case, namely from Eqs. (102) and (103) respectively. 

Figure Q shows profiles of the density perturbation <j(^i^ and the momentum density perturbation Ji for the close- 
limit, boosted, head-on collision. Recall that for inspiral and head-on collisions Ji = 2 J2 and J3 = 0. The perturba- 
tions (7(1) and Ji were calculated from neutron stars with initial separation ^ = 0.1i?(o) and velocity v = 0.1c. The 
corresponding perturbations for the inspiral case only differ from t he p erturbations shown in Fig. ^ by a constant 
numerical factor (compare Eqs. (|9^ ) and (101) with Eqs. (109) and (IH). The TOV parameters for the background 

and colliding stars are: = 2.69 x lO^^g/cm"^ and K(o) = 100 km^. For these parameters, the mass and radius of the 
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background star are M(o) = 1.24 Mq and i?(o) = 9.0 km, respectively. The initial colliding stars, which are displaced 

a distance 0.1i?o from the center of mass, follow from p* = 2.98 x lO^^g/cm^ and = 90.25 km^ With these 
parameters, the colliding stars have a mass and radius of Af* = 1.17 and i?o* = 8.58km, respectively. Figure || 
shows the solutions to the conformal perturbation and the harmonic components Wi and W2 (see Eq. ^) of the 
vector for the close-limit collision of neutron stars corresponding to the perturbations in Fig.nl 



VI. CONCLUDING REMARKS 



In this paper, we have presented a framework for constructing initial data relevant for perturbative studies of neutron 
stars. Our approach was to "linearize" Lichnerowicz- York's standard procedure for the initial- value problem in General 
Relativity, and it facilitates (to a certain extent) setting astrophysical initial data for perturbation evolutions, cf. 
It is straightforward to compare our method (as well as the results) to the fully nonlinear one, which is important 
since a main motivation for perturbation studies is to provide benchmark tests for nonlinear numerical relativity. 

As examples of interesting initial data that can be constructed from our equations, we constructed data for merging 
neutron stars in the close-limit approximation. The simplest case of these data sets, that describe head-on collision of 
two initially static stars, has already been extensively discussed in No studies of the more general data with initial 
momentum and for inspiralling collisions have yet been performed. Such simulations should obviously be carried out, 
and we hope to be able to discuss the relevant results, as well as possible extensions of the framework developed in 
this paper to, for example, rotating configurations, in the near future. 
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FIG. 1. Density and momentum density Ji perturbations (see Eqs. and 



101 in the text) for the close-hmit, boosted, 
head-on colhsion of neutron stars with initial separation ^ = 0.1 7?(o) and velocity v = 0.1 c. The corresponding perturbations 
for the inspiral case differ from these quantities by constant numerical factors. The TOV parameters for the background and 
colliding stars are: pc"' = 2.69 x 10^^ g/cm"^ and K(o) = 100 km^. For these parameters, the mass and radius of the background 
star are M(o) = 1.24 Mq and i?(o) = 9.0 km, respectively. The initial colliding stars, which are displaced a distance O.liio from 
the center of mass, follow from p* = 2.98 x 10^^ g/cm^ and — 90.25 km^. With these parameters, the colliding stars have a 
mass and radius of Af, = 1.17 M© and Ro* = 8.58km, respectively. 
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